BIHEP-EP1_98-01 
UH-51 1-903-98 

Branching fractions for ifj(2S) — ► 777' and 777 



J. Z. Bai 1 , J. G. Bian 1 , I. Blum 11 , Z. W. Chai 1 , G. P. Chen 1 , H. F. Chen 10 , J. Chen 3 , 
J. C. Chen 1 , Y. Chen 1 , Y. B. Chen 1 , Y. Q. Chen 1 , B. S. Cheng 1 X. Z. Cui 1 , H. L. Ding 1 , 
L. Y. Ding 1 , L. Y. Dong 1 , Z. Z. Du 1 , W. Dunwoodie 7 , S. Feng 1 , C. S. Gao 1 , M. L. Gao 1 , 

S. Q. Gao 1 , P. Gratton 11 , J. H. Gu 1 , S. D. Gu 1 , W. X. Gu 1 , Y. F. Gu 1 , Y. N. Guo 1 , 
S. W. Han 1 , Y. Han 1 , F. A. Harris 8 , J. He 1 , J. T. He 1 , M. He 5 , D. G. Hitlin 2 , G. Y. Hu 1 , 
H. M. Hu 1 , J. L. Hu 1 , Q. H. Hu 1 , T. Hu 1 , X. Q. Hu 1 , J. D. Huang 1 , Y. Z. Huang 1 , 
J. M. Izen 11 , C. H. Jiang 1 , Y. Jin 1 , Z. J. Ke 1 , M. H. Kelsey 2 , B. K. Kim 11 , D. Kong 8 , 
Y. F. Lai 1 , P. F. Lang 1 , A. Lankford 9 , C. G. Li 1 , D. Li 1 , H. B. Li 1 , J. Li 1 , P. Q. Li 1 , 
R. B. Li 1 , W. Li 1 , W. D. Li 1 , W. G. Li 1 , X. H. Li 1 , X. N. Li 1 , H. M. Liu 1 , J. Liu 1 , 
J. H. Liu 1 , R. G. Liu 1 , Y. Liu 1 , X. C. Lou 11 , B. Lowery 11 , F. Lu 1 , J. G. Lu 1 , J. Y. Lu 1 , 
L. C. Lu 1 , C. H. Luo 1 , A. M. Ma 1 , E. C. Ma 1 , J. M. Ma 1 , R. Malchow 3 , H. S. Mao 1 , 
Z. P. Mao 1 , X. C. Meng 1 , J. Nie 1 , S. L. Olsen 8 , J. Oyang 2 , D. Paluselli 8 , L. J. Pan 8 , 
J. Panetta 2 , F. Porter 2 , N. D. Qi 1 , X. R. Qi 1 , C. D. Qian 6 , J. F. Qiu 1 , Y. H. Qu 1 , 
Y. K. Que 1 , G. Rong 1 , M. Schernau 9 , Y. Y. Shao 1 , B. W. Shen 1 , D. L. Shen 1 , H. Shen 1 , 
X. Y. Shen 1 , H. Y. Sheng 1 , H. Z. Shi 1 , X. F. Song 1 , J. Standifird 11 , F. Sun 1 , H. S. Sun 1 , 
S. Q. Tang 1 , W. Toki 3 , G. L. Tong 1 , F. Wang 1 , L. S. Wang 1 , L. Z. Wang 1 , M. Wang 1 , 
Meng Wang 1 , P. Wang 1 , P. L. Wang 1 , S. M. Wang 1 , T. J. Wang 1 !, Y. Y. Wang 1 , 
M. Weaver 2 , C. L. Wei 1 , Y. G. Wu 1 , D. M. Xi\ X. M. Xia 1 , P. P. Xie\ Y. Xie 1 , 
Y. H. Xie 1 , W. J. Xiong 1 , C. C. Xu 1 , G. F. Xu\ S. T. Xue 1 , J. Yan 1 , W. G. Yan 1 , 
C. M. Yang 1 , C. Y. Yang 1 , J. Yang 1 , W. Yang 3 , X. F. Yang 1 , M. H. Ye 1 , S. W. Ye 10 , 
Y. X. Ye 10 , K. Yi \ C. S. Yu 1 , C. X. Yu 1 , Y. H. Yu 4 , Z. Q. Yu 1 , Z. T. Yu 1 , C. Z. Yuan 1 , 
Y. Yuan 1 , B. Y. Zhang 1 , C. C. Zhang 1 , D. H. Zhang 1 , Dehong Zhang 1 , H. L. Zhang 1 , 
J. Zhang 1 , J. L. Zhang 1 , J. W. Zhang 1 , L. S. Zhang 1 , Q. J. Zhang 1 S. Q. Zhang 1 , 
X. Y. Zhang 5 , Y. Zhang 1 , Y. Y. Zhang 1 , D. X. Zhao 1 , H. W. Zhao 1 !, J. W. Zhao 1 , 
M. Zhao 1 , W. R. Zhao 1 , Z. G. Zhao 1 , J. P. Zheng 1 , L. S. Zheng 1 , Z. P. Zheng 1 , 
G. P. Zhou 1 , H. S. Zhou 1 , L. Zhou 1 , Q. M. Zhu 1 , Y. C. Zhu 1 , Y. S. Zhu 1 , B. A. Zhuang 1 

(BES Collaboration) 
1 Institute of High Energy Physics, Beijing 100039, People's Republic of China 
2 California Institute of Technology, Pasadena, California 91125 
^Colorado State University, Fort Collins, Colorado 80523 
4 Hangzhou University, Hangzhou 310028, People's Republic of China 
5 Shandong University, Jinan 250100, People's Republic of China 
6 Shanghai Jiaotong University, Shanghai 200030, People's Republic of China 
7 Stanford Linear Accelerator Center, Stanford, California 94309 
8 University of Hawaii, Honolulu, Hawaii 96822 
9 University of California at Irvine, Irvine, California 92111 
10 University of Science and Technology of China, Hefei 230026, People's Republic of China 
11 University of Texas at Dallas, Richardson, Texas 15083-0688 



1 



Abstract 



We report first measurements of the branching fractions: B(ifj2S ~^ iv') = 
(1.54±0.31±0.20) x 10" 4 and B{ip 2S -» 7v) = (0.53±0.31±0.08) x 10~ 4 . The 
ip(2S) — > 77?' result is consistent with expectations of a model that considers 
the possibility of rj-r\ c mixing. The ratio of the ip(2S) — > 77/ and ip(2S) — ► 777 
rates is used to determine the pseudoscalar octet-singlet mixing angle. 

I. INTRODUCTION 

In contrast to the J/ip, experimental results on radiative decays of the if)(2S) to non- 
charmonium hadrons are scarce; in the latest Particle Data Group tables, only upper limits 
for a few decay modes are listed [|TJ. Moreover, it is found experimentally that, while decays 
to prr and K*K Vector Pseudoscalar (VP) final states are significant (~ 1%) for the J/ip, 
hadronic decays of the ip(2S) to these same VP final states are strongly suppressed P,H. 
This longstanding mystery of charmonium physics is referred to in the literature as the pir 
puzzle M. The processes J/if> — > 777' (958) and 777 are radiative VP channels that have been 
measured by several experiments 0. It is of interest to see if the same radiative VP decays 
of the ip(2S) are suppressed to the same extent as the hadronic pit and K*K decays. 

Recently, the CLEO experiment has reported an anomalously large branching fraction 
for the inclusive production of rf in the 5-meson decay B — > t]'X s , where X s denotes an 
inclusive hadronic system containing a strange quark H . One possible interpretation is the 
presence of an intrinsic charm component of the 77' meson induced by the strong coupling 
of the 7/ to gluons via the QCD axial anomaly |5]]. The resulting r(-r\ c mixing has been 
proposed as the dominant mechanism for the OZI forbidden radiative charmonium decays 
such as tp(nS) — ► 777' and 777. In this case, the branching fraction for ip(2S) — > 71]' is 
estimated to be in the range (1.0 — 2.7) x 1CT 4 0. 

The ratio of the 777' and 777 decay rates of the J PC = 1 charmonium states is sensitive 
to the pseudoscalar octet-singlet mixing angle 9 p . Assuming that the process occurs primarily 
through radiation of the photon from one of the initial state c-quarks and the applicability 
of SU(3) symmetry for the decay amplitudes, one has the simple relation H 

r(^/y? -» ij) _ fprA 3 1 

r(Vv? -> iv) \PvJ tan 2 6» p ' " 

where p v (jy) is the momentum of the rj (?/) in the ip(nS) rest frame. The measured J/ip 
branching fraction values [J] imply a mixing angle of \9 P \ = 22° ± 1° ± 4°, (the second 
error is theoretical), which agrees well with the value determined from other processes ||. 
Measurements of the corresponding branching fractions for the ip{2S) provide a consistency 
check of this relation. 

In this report we present the first measurement of the branching fractions for ip(2S) — > 777' 
and ip(2S) — > 777 using 3.7x 10 6 ip{2S) decays collected using the Beijing Spectrometer (BES) 
located at the Beijing Electron-Positron Collider (BEPC) at the Beijing Institute of High 
Energy Physics. 
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II. THE BES DETECTOR 

The BES is a large solid-angle magnetic spectrometer that is described in detail in ref . |TD[ . 

Charged particle momenta are determined with a resolution of a p /p = 1.7%yl + p 2 (GeV 2 ) 
in a 40-layer cylindrical drift chamber. Radially outside of the drift chamber is a 12-radiation- 
length barrel shower counter (BSC) comprised of gas proportional tubes interleaved with 
lead sheets. The BSC measures the energies and directions of photons with resolutions of 
oeIE ~ 22%/ J 1 E(GeV), = 4.5 mrad, and a$ = 12 mrad. The iron flux return of the 
magnet is instrumented with three double layers of counters that are used to identify muons. 

For this analysis we use charged tracks with momentum greater than 80 MeV/c that 
are well fit to a helix originating near the interaction point. Candidate 7's are associated 
with energy clusters in the BSC that have more than three hit tubes in at least two readout 
layers. We use charged tracks and 7's that are within the polar angle region | cos#| < 0.8. 
We reject tracks that are identified as muons, or that produce high energy showers in the 
BSC that are characteristic of electrons. When computing energies, each charged track is 
assigned the pion mass. 



III. DETERMINATION OF B(tp(2S) -> 7/7') 

For the ip(2S) — > 777' measurement we investigate the decay chains 
^p2S ->• 7 V' and 1P2S -> 7 V' 

7T + 7T 77 

u 

7T + 7T U27. 

It follows that the reactions of interest are ip(2S) — ► 7r + 7r~77 for the p°7 mode, and ip(2S) 
7r + 7r~777 for the 7r + 7r _ r/ decays. 



A. The ip(2S) — ► 7/7' — > 77P° measurement 

For the measurement using the rj' —>■ 7p° mode, we require two oppositely charged tracks 
with an opening angle 8 open < 130° and at least two candidate 7's that are more than 10 
degrees away from the nearer charged track. The events where the total energy of the two 
charged tracks is less than 2.1 GeV are subjected to a four-constraint kinematic fit to the 
hypothesis ip(2S) — > 7r + 7r~77, and required to have x 2 < 15. The 7r + 7r~7 mass distribution 
for events with M n + n - within 0.15 GeV of M p and a 77 opening angle greater than 110° is 
plotted in Fig. [l], where a peak at the mass of the rj' (958) is apparent. 

The curve in Fig. [1] is the result of a fit to the measured mass distribution with the 
r;(958) represented as a Gaussian, and a third-order polynomial background function. The 
width of the Gaussian is fixed at the Monte-Carlo determined experimental resolution of 



a = 0.01 GeV fill . The fitted Gaussian has a peak position at M f? ( 958 ) and N evts = 28.1 ±7.2 
events. Events from the cascade decays ip(2S) — ► anything + J/ip, where J/ip — > 77/ or 
J Jijj — > vr°p°, also can give a peak at M^c^y We subject a sample of Monte Carlo simulated 
events equivalent to ten times the i[>(2S) data set to the same selection and fitting procedure. 
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FIG. 1. The 7r + 7r 7 invariant mass distribution for selected events. 



The resulting estimate of the contamination from this source is N^ g = 1.4±0.5 events, where 
the error is statistical and comes from the fit. 

We use Monte Carlo simulated events to determine the acceptance. The events are 
generated with a 1 + cos 2 6 angular distribution for the ip(2S) — > 777' decays, an isotropic 
distribution for the rj' —>■ 7p° decays, followed by helicity ±1 p° — > 7r + 7r~ decays. The 
acceptance determined in this way is e P7 = 0.18 ± 0.02, where the error includes both 
Monte Carlo statistics (7%) and uncertainties in the simulation program (8%) 
quadrature. 

The '0(25') — * 777' branching fraction is determined from the relation 



121 added in 



2S 



N, 



evts 



N, 



bkg 



(1.36 ±0.37 ±0.20) x KT 4 . 



(2) 



Here the first error is statistical, and the second is the systematic error due to uncertainties 
in Nfag (9%), the acceptance (11%) and the 7/ — > 7p branching fraction (3.3%) added in 
quadrature. 



B. The tjj(2S) — ► 7/7' — *■ 777+77 77 measurement 

For the measurement using the rf — > 7r + 7r _ 77 mode, we require two oppositely charged tracks 
with an opening angle open < 70°, and at least three candidate 7's that are more than five 
degrees away from the nearer charged track. We select events with a total energy for the 
two charged tracks that is less than 1.2 GeV, and require them to satisfy a four-constraint 
kinematic fit to the hypothesis ip{2S) — > 7r + 7r~777 with x 2 < 12. We identify 77 pairs with 
an invariant mass within 0.03 GeV of M v as candidate 77's. The 7i + 7i^i] mass distribution 
for the selected events is plotted in Fig. ^|. There is a peak in the data at the mass of the 
7/(958). 

The curve in Fig. |2] is the result of a fit to the measured mass distribution with the 7/(958) 
represented as a Gaussian, and a polynomial background function that is forced to zero at 
the 7r + 7r~77 threshold. The width of the Gaussian is fixed at a = 0.018 GeV, the resolution 
value determined from the MC simulation. The fitted Gaussian has N evts = 16.8 ± 4.9 
events. The Monte Carlo estimate of backgrounds from ip(2S) —>■ J/ip cascade decays is 
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FIG. 2. The Tr + ir~rj invariant mass distribution for selected events. 

Nbkg = 0.35 ± 0.03. The MC determined acceptance for this mode is 0.14 ± 0.015, and the 
corresponding ip(2S) — > 77/ branching fraction is 

B{ip 2 s -> iv) = (2-00 ± 0.59 ± 0.29) x 10~ 4 . (3) 
The agreement with the result determined for the 7/ — > 7p° mode is reasonable. 

IV. THE ip(2S) -s- 777 MEASUREMENT 

For the ip{2S) — > 717 measurement, we use the 77 — > tt + 7t~7t° decay mode. This corresponds 
to the same ip{2S) — > 7r + 7r _ 777 reaction as for the 77' — > 7T + 7T~i] measurement. 

We require two oppositely charged tracks with total energy less than 1.7 GeV, and at 
least three candidate 7's that are more than five degrees away from the nearer charged 
track. The events are required to satisfy a four-constraint kinematic fit to the hypothesis 
ip(2S) — ► 7r + /T~777 with % 2 < 12. We identify 77 pairs with an invariant mass within 
0.025 GeV of m^o as candidate vt°'s. The 7r + /T~7r° mass distribution is plotted in Fig. |3|. 
Here small clusters of events appear at the mass of the ?7(547) and the u(780). 

6 - 




M(n+n-nO)' (GeV) 

FIG. 3. The ir + iT~iT invariant mass distribution for selected events. The region near the 
cj(785) is excluded from the fit. 

The curve in Fig. || is the result of a fit to the measured mass distribution with the 
7] represented as a Gaussian, and a polynomial background function. The width of the 
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Gaussian is fixed at a = 0.013 GeV, the value determined from the MC simulation. The 
a;(780) mass region is excluded from the fit. The fitted Gaussian has a peak position at 
M^+^-^o = 0.56 ± 0.01, which is one standard deviation above M n , and an area of N^ vts = 
4.1 ± 2.4 events. 

As a check, we used the events in the 7r° sidebands of the 77 invariant mass distribution 
as an experimental estimate of our background. Here we find no events within ±3a of M v 
and a fit to the sideband-subtracted M n+7r _ n o distribution yields 6.0 ± 2.5 r/ events. 

The ip{2S) — > 77/ signal has a statistical significance corresponding to a little less than 
2a P"3|| . If we treat the 4.1 observed events as a real signal, the if>(2S) — > 777 branching 
fraction is determined to be 

Bfas -> m) = (0.53 ± 0.31 ± 0.08) x 10" 4 . (4) 

(The MC-determined acceptance for this channel is 0.10 ± 0.012.) The 4.1 events from the 
fit imply a 90% confidence level (c.l.) upper limit of 7.2 events; this corresponds to a 90% 
c.l. limit on the ip(2S) — > 777 branching fraction of 0.9 x 10 -4 . 



V. DISCUSSION 

Combining the two results for B(ip(2S) — > 777') from the different 7/ decay modes gives [|TJ] 

B{ip 2 s -> li) = (1-54 ± 0.31 ± 0.20) x 10" 4 , (5) 

which is within the range expected for the case where 77' — r\ c mixing is important 0. To 
compare with J /if} decays, we use the ratio 

Qvt = WTF, ~ = 0-036 ± 0.009. (6) 

This low value for Q 7r) / indicates that this ip(2S) decay mode is suppressed relative to di- 
lepton decays, where the corresponding ratio Qu = 0.147 ±0.023 but not as severely as 
in the case of pit, where Q p7T < 0.002, or K* + K~, where Q k *+k- < 0.006 @. Pinsky |I5| 



relates the processes ip(2S) — > ^vj to the hindered Ml transition ip(2S) — ► 7/7 c . He predicts 
Q 1V > = 0.002, which is well below our measured value. 

The suppression of J/tf> — > 777 relative to J /if) 777' decays appears to also occur for 
the ip(2S): 

B(J/if> - 1V ) 



B(J/if> <yrf) 
B(iP(2S) - 777) 



0.200 ± 0.023 (PDG); (7) 
0.34 ±0.22 (this expt). (8) 



B(ip(2S) -> 777' 

In nrinn /~\4- 4-^-1^ 1 1 1 t Tn -Ki ^-r r\ v-i rr Ia at I i , . , _ 

y p | — z,o _ 10 o, 



Our results provide an independent evaluation of the mixing ang le of |0J = 28°±Ioo, whi ch 



is consistent with other determinations, albeit with larger errors. 
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